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Abstract 

We outline two cubic spline based approximation schemes for the 
estimation of structural parameters associated with the transverse vibration 
of flexible beams with tip appendages. The identification problem is 
formulated as a least squares fit to data subject to the system dynamics which 
are given by a hybrid system of coupled ordinary and partial differential 
equations. The first approximation scheme is based upon an abstract semigroup 
formulation of the state equation while a weak/variational form is the basis 
for the second. Cubic spline based subspaces together with a Rayleigh-Ritz- 
Galerkin approach was used to construct sequences of easily solved finite 
dimensional approximating identification problems. Convergence results are 
briefly discussed and a numerical example demonstrating the feasibility of the 
schemes and exhibiting their relative performance for purposes of comparison 
is provided. 
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In this short paper we briefly outline two cubic spline based 
approximation schemes for the solution of inverse problems involving the 
vibration of flexible beams with attached tip bodies. The identification 
problem is formulated as the least squares fit to data of a hybrid system of 
coupled partial and ordinary differential equations describing the dynamics of 
the beam and tip bodies. The resulting optimization problem is infinite 
dimensional and as such, necessitates the use of some form of approximation. 
The schemes we have developed are based upon the construction of a sequence of 
approximating identification problems in which the underlying constraining 
state equations are semi-discrete finite dimensional approximations to the 
infinite dimensional distributed system which governs the original 
identification problem. Our study includes both theoretical convergence 
results and numerical testing. 



Although our general approach applies to a broad class of problems, to 
illustrate the two methods, we consider the problem of identifying the 
spatially invariant flexural stiffness q x and linear mass density q 2 for a 
beam of length l clamped at one end and cantilevered at the other with a tip 
(point) mass of magnitude q3 which is also to be identified. Using the 
Euler-Bernoulli theory and elementary Newtonian mechanics, the system of 
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q 2 u tt (t ' X> * - q l"xxxx (t > x) + f(tlX) 
«3 u tt (t>t) ' q l u xxx <t ’ ll> + 8<t) 


X e (0 ,Jt) , t e (0,T) (1) 

t e (0,T) (2) 


describing the transverse deflection of the beam and tip mass is obtained 
where f and g denote externally applied lateral loads (see [2]). The boundary 
conditions are given by 


u(t,0) = u x (t,0) = u xx (t,4) = 0 


t e (0 ,T) , 


(3) 


The initial conditions are of the form 


u(0,x) = u t (0,x) = i|»(x) 


x e [0 ,£] 


(4) 


The identification problem is stated formally as 


T +3 

(ID): Find q = (q 1 ,q 2 »q 3 ) e Q, Q a compact subset of R , 

which minimizes 

j(q;u(*,«;q)) = ]£ £ |u(t ,x ;q) - u(t ,x )| 2 
i=l j=l J J 

subject to u(»,»;q) 

{u(t ± *Xj ) } i=i , . . v 
j = l * • • • > V 

observations. 


being a solution to (1) - (4) where 
denote a set of given displacement 



The first approximation scheme, which is described in detail in [6], is 
based upon the recasting of (1) - (4) as an abstract evolution equation set in 
an infinite dimensional Hilbert space. Let H = RxH°(0,fc) with inner product 
<(n,<f>) ,(C ,t)>H “ HC + where {h^ c ,<*,*>j c } denote the usual Sobolev 

spaces and Sobolev inner products. Define the operators MQ(q):H-*-H and 
A 0 (q):DCH*H by M 0 (q)(n,<()) = (q 3 n,q 2 <|>) and 

,<fr) = C — q^4» (i) respectively where 

D = {(n,<J>)eH:<freH 4 (0,£), <(>(0) = <T(0) - 0, 4>"(l) = 0, q = 4>U)}. Assume 
f e L 2 ([0,T],H°(0,i)), g e L 2 (0,T), $ e H 2 (0,£) and e H°(0,£) and let 
F 0 (O “ (g(t)»f(t,*)) , $ = (<K*),<I>) and i|> = ('|>(£),'j0 where ip(i) e R is 

f 

specified if it is not well defined. We then rewrite (1) - (4) as an abstract 
second-order system in H; 

MgCq) D^u(t) + AQ(q)G(t) = F Q (t) , t e (0,T) (5) 

u(0) = l, D t G(0) = * (6) 

where u(t) = (u(t,£) ,u(t,-)) e H and D t - ar • 

V - {(n,«.)eH:$eH 2 (0, »),$(<>) = ♦'«>) . 0,n - ««) and define 

L:VCH+H by L( ,+) = (0,$"). The operator Ag(q) can then be written in 
factored form as AQ(q) = qjL*L where L* :Dom(L*) C H+H is given by 
L*(n,<f>) = (-<)>'(£), <f>") for (n,<f>)eDom(L*) = {(n,<|>)eH : <|>eH 2 (0,£), <f.< «.) = o}. 

Let Z = H x H with inner product 

<z>w>q = ^ q l z i» w i' > H + Choosing our state as 

z(t) = (Lu(t) ,D t u(t)) we rewrite (5), (6) as the first-order system in Z given 


by 
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z(t) = A(q)z(t) + F(t;q), z(0) = z Q 


(7) 


where A(q) :Dom(L*)xVCZ-»-Z is given in matrix form by 


A(q) = 


0 

-q^Cq) x l* 


L 

0 


F(t;q) = (0,M 0 (q) _1 F 0 (t)) and zq - (L$,$) ($ is assumed to be an element in 

V). The operator A(q) is densely defined and skew self adjoint. It 

therefore generates (see [8]) a Cq group of unitary operators 

(S(t;q):^» < t < »} on Z. The mild solution (strong or classical if F and 

Zq are sufficiently regular) to (7) is then given by (see [5]) 


t 

z(t) = S(t;q)z_ + / S(t-x;q)F(x;q)dx. (8) 

U 0 


Problem (ID) can then be written as 


(IDA): Find qeQ which minimizes j(q;Ci( «)z( »;q)) subject to 

z(*;q) being given by (8) where the operators 
C^(x):Z-»-R are defined by 

x x 

c, (x)((n, ♦),(?,♦)) = 11 <K<j)dodx for x e [0,*]. 

1 0 0 

N 

For each N = 1,2..., let S^(A ) denote the space of cubic spline 
functions on the interval [0,£] corresponding to the uniform partition 
A N ={°4, (see [7]). Let W N = { ( n , <f>) eH : «J>eS 3 ( A N ) , «K A) = 0 } , 

V N = {(n,4>)eH:<M:S 3 (A N ),<|>(0) = <T(0) - 0,n = ♦(!)} and Z N = W N x V N . 


Then 
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Z N is a 2N+4 dimensional subspace of Z which satisfies 

Z^CDom(A(q)) = Dom(L*)xV for all qeQ. Let denote the orthogonal 
projection of Z onto Z^ with respect to the <•,•> inner product. 
Define A^(q):Z^ -*• Z* 1 by A^(q) = P^A(q) and noting that A^(q) is a linear 
operator defined on a finite dimensional space (and is therefore bounded), 
consider the initial value problem 


z N (t) = A N (q)z N (t) + P^(t)F(t;q), t £ (0,T) 


z N (0) = P*z Q , 


and its solution 

N, v f A N, x a N 

z (t) = exp( A (q)t)P^z Q 

+/ exp[A N (q)(t-T))p^F(T;q)dx t e [0,T]. (9) 

0 q 

The approximating identification problems for the first method take the form 

(IDN!): Find qeQ which minimizes j(q;C^( *)z N (*;q)) subject 

to z N (»;q) being given by (9). 

Under rather mild assumptions it can be shown that each of the problems 

(IDNj) admits a solution q^ and that the sequence {q^} contains a 

Nu _N k _ 

convergent subsequence lq }, q + q e Q as k + ®. Using standard 

approximation results from linear semigroup theory (see [5]) and the 
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properties of spline functions (see [7]) to establish convergence of the state 

approximations, it can then be argued that cf is a solution to problem (ID). 

The second approximation scheme involves the rewriting of (1) - (4) in an 

equivalent weak/ variational form (see [4]). Let H and V be as they have 

been defined previously. Define the inner product on V, < # , # >y by 

= <L<}> ,L$>„. We have the usual dense embeddings V C H C V' where V" 
V H 

is the dual of V. We consider the weak form of (1) - (4) 

<M o (q)D^u(t),0> H + aCu(t),0;q) = <F Q (t),0> H> 0 e V, t e (0,T) (10) 

u(0) = <f>, D t u(0) = (11) 

where the bilinear form a(*,*;q) on VxV is defined by 

A A A A A A 

a(<fr,ij>;q) = q^ <L<j> ,L iJj> h = <AQ(q)<J>,i{>> H . The derivatives in the definition of 
the operator Ag(q) are interpreted in the distributional sense and Ag(q) 
is considered to be an element in Jz^(V,V'), the space of continuous linear 
operators from V into V'. 

We define a Galerkin approximation u^ to u as the solution to 

<M o (q)D^u N (t),0 N > H + a(G N (t),0 N ;q) = <F Q (t),0 N > H , 0 N e V N , t e (0,T) (12) 

u N (0) = P N $, D t u N (0) = P N i (13) 

where is as it was defined above and P^ is now the orthogonal 

projection of H onto with respect to the i nner product. 
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The approximating identification problems are then given by 

(IDN 2 ): Find qeQ which minimizes j(q;C 2 ( «)u^( »;q)) subject to 

u (*;q) being the solution to (12), (13) where the operators 
C 2 (x):V+R are defined by C 2 (x) (<|>(A) ,$) = <|»(x) ,xe[0, A] . 

Using standard variational arguments of the type found in [3] it can be 
argued that under sufficient regularity assumptions, for (q^}CQ with 
qN+q* we h ave u N (q N )-»-u(q*) in V and D t u N (q N )+D t u(q*) in H as 

N-*» where u^(q^) is the solution to (12), (13) corresponding to q^ and 
u(q*) is the solution to (10), (11) corresponding to q*. This in turn 

yields a convergence result analogous to the one stated above for the first 
scheme. A more detailed discussion of these results appears in [1]. 

We demonstrate the feasibility of our methods with a simple example. We 
consider a beam of length 1 and use the two schemes described above to 
estimate its stiffness qj, its linear mass density q 2 and the magnitude of 
a tip mass qg. We assumed that the system was initially at rest (<£=i(f= 0) 
and then excited via the distributed lateral load f(t,x) = e x sln 2 irt and the 

—9 . 

point force g(t) = 2e“ acting on the tip mass. Displacement observations 
at positions Xj = .5, .75, 1.0, j=l,2,3, at times t ± = .5i, i=l ,2, . . . , 10, 
were generated using the "true" parameter values qj = 1 . 0 , q 2 = 3.0 and 
<13 = 1«5, the first two natural modes of the system and a standard Galerkin 
method. The approximating identification problems (IDN^) and (IDN 2 ) were 
solved using an iterative Levenberg-Marquardt scheme with the approximating 
state equations solved at each iteration using a variable order Adams 
Predictor Corrector method. The "start-up" values for the unknown parameters 
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were taken to be q® = .7, = 2.7 and q^ = 1.7. Our results for methods 1 

and 2 are given in Tables 1 and 2 respectively below. 


Table 1 


N 

-N 

q l 

-N 

q 2 

-N 

q 3 

CPU (m: s) 

2 

1.00083 

2.99817 

1.50090 

0:08.70 

3 

1.00072 

3.00317 

1.49868 

0:29.38 

4 

1.00061 

2.99227 

1.50141 

1:01.19 

5 

1.00009 

2.98711 

1.50195 

1:35.26 

6 

1.00061 

2.99039 

1.50144 

2:45.20 


Table 2 

N 

-N 

q l 

-N 

q 2 

-N 

q 3 

CPU (m: s) 

2 

1.00057 

3.04455 

1.48957 

0:09.19 

3 

1.00067 

3.01256 

1.49707 

0:22.10 

4 

1.00027 

3.00922 

1.49721 

0:57.57 

5 

1.00016 

2.98936 

1.50262 

1:22.52 

6 

.99912 

2.99720 

1.49952 

2:52.76 
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